Fast electrochemical doping due to front instability in organic semiconductors 
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The electrochemical doping transformation in organic semiconductor devices is studied in applica- 
tion to light-emitting cells. It is shown that the device performance can be significantly improved by 
utilizing new fundamental properties of the doping process. We obtain an instability, which distorts 
the doping fronts and increases the doping rate considerably. We explain the physical mechanism of 
the instability, develop theory, provide experimental evidence, and perform numerical simulations. 
We further show how improved device design can amplify the instability thus leading to a much 
faster doping process and device kinetics. 



Organic semiconductors (OSCs) are expected to revo- 
lutionize everyday electronics by offering interesting and 
attractive properties which distinguish them from con- 
ventional inorganic semiconductors [H-lsj. In addition to 
simple processing, low cost, soft and conformable char- 
acter, OSCs provide the intriguing possibility of in-situ 
chemical and electrochemical doping. This doping trans- 
formation leads to significant changes in important ma- 
terial properties, such as conductivity, color, volume, and 
surface energy [l, Q . The opportunity for a controllable 
tuning of the properties of OSCs via doping has stim- 
ulated the emergence "organic electronics" with a large 
number of novel and flexible applications Opera- 
tionally, electrochemical doping is performed by apply- 
ing a potential to a metal electrode coated with an OSC 
in contact with an electrolyte. During n-type (p-type) 
doping, electrons (holes) are injected into the OSC from 
the cathode (anode) and subsequently compensated by 
redistribution of cations (anions) in the electrolyte [5|-l9| . 
When OSC is populated by light charges, its conductivity 
increases by two-three orders of magnitude in comparison 
with the original conductivity due to bulky ions. A com- 
plex manifestation of the electrochemical doping is the 
in-situ formation of a dynamic p-n junction when p- and 
n-doping fronts meet in a light-emitting electrochemical 
cell (LEC) ^^E^. However, the electrochemical dop- 
ing involves redistribution of heavy and slow ions, which 
makes such organic electronics devices slow to turn-on. 
The speed-up of the doping process is thus vitally im- 
portant for technical applications. Here we obtain, the- 
oretically and experimentally, instability of the doping 
transformation front, which leads to much faster kinetics 
of the process. We suggest an improved device design on 
the basis of the attained fundamental insight. We show 
that a corrugated pattern in the electrode surface am- 
plifies the instability and increases the doping rate; thus 
substantially reducing the device turn-on time. 

The configuration of two doping fronts counter- 
propagating towards each other in an LEC device is 
schematically illustrated in Fig. 1(a). Fig. 1(b) shows 



the internal structure of two such planar fronts, as cal- 
culated using the method of [19]. The doping fronts, 
however, can never be ideally planar because the unavoid- 
able small perturbations grow with time and distort their 
shape. The physics behind this new instability is related 
to the phenomenon of St. Elmo's fire, caused by the in- 
creased electric field at a convex conducting surface. In 
the present configuration, any leading perturbation hump 
at a doping front causes local increase of the electric field 
in the undoped region. Since the front velocity is pro- 
portional to the electric field, see Eq. (1) below, this 
stronger electric field will force the front to propagate 
faster locally, thus producing a positive feedback and an 
unstable growth of the hump. 

To study the instability growth, we employ the p/n- 
doping front velocity [HI 
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where Eq = — V^o is the electric field intensity just ahead 
of the front, the electric field potential, no the ini- 
tial ion concentration in the pristine active material, nh,e 
the final concentration of the holes/electrons behind the 
front (determined by the thermodynamic properties of a 
particular OSC), /i+ and /i_ the mobilities of the posi- 
tive/negative ions. The " ±" sign in Eq. (1) indicates that 
the p-type (+) and n-type (— ) doping fronts propagate in 
opposite directions. We consider a perturbation Fourier 
mode X{t) ex.p{iky) of a p-doping front x = Xp{y^t) in 
Fig. 1(a), which induces also perturbations of the elec- 
trical field in the undoped region. The doped region may 
be treated as equipotential due to high conductivity. The 
front thickness may be characterized by a length scale Lp, 
of about 10~^ — 10~^mm, determined by ionic diffusion 
[Toj . Since the characteristic size of the experimentally 
observed front perturbations (approximately 0.2mm) is 
much larger than Lp, then we may treat the front as in- 
finitesimally thin. Most of the time the doping fronts are 
sufficiently far away from each other, k{Xn — Xp) 1, 
therefore, the instability of one front is not affected by 
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FIG. 1. Schematic of the p- and n-doping fronts in an organic 
semiconductor film (a), and internal structure of the planar 
fronts at 90 s (b). 



the other. The linearized form of Eq. (1) 
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together with the solution to the Laplace equation for the 
electric potential in the udoped region, (j) ex ex.p{iky — 
\k\x), yields the equation dtX = Up\k\X describing the 
perturbation growth with time. Thus, we obtain an in- 
stability for the p-type doping front, where the initially 
small perturbations grow exponentially as X ex exp(crt), 
with the growth rate a = Up\k\. The same result holds for 
the n-type doping front by replacing Up with Un- This 
new instability shows interesting mathematical similar- 
ities to the Darrieus-Landau instability encountered in 
combustion [2o|, Hlj and inertial confined fusion [22]. 

To observe the instability experimentally, we utilize an 
open planar LEG device, comprising an {MEH-PPV + 
PEG + KGF3SG3} active material positioned between 
two Au electrodes as described in [19|. Two counter- 
propagating doping fronts can be distinguished in Fig. 
2 (a, b); the dark regions with quenched fluorescence 
correspond to the doped MEH-PPV; the electrodes are 
indicated with white dashed lines. When the fronts meet 
and form a p-n junction, recombination of subsequently 
injected holes and electrons can lead to light emission 
(see Fig. 2(c), taken in darkness). Figures 2 (a,b) clearly 
demonstrate that both doping fronts are unstable with 
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FIG. 2. Experimental photos of the doping fronts demon- 
strate development of the instability at the initial stage (a), 
t = 46 s, at the developed stage (b), t = 80 s, and the dy- 
namic p-n junction (c), t > 170 s. Plot (d) compares the 
p-front shapes obtained in the simulations (white curve) and 
experiments (shading) at 70 s. Plots (e), (f) show relative 
increase of the electric field in the undoped region obtained 
numerically for the whole front squeezed along Y-axis, (e), 
and for the selected part with equal scales, (f). 



respect to small perturbations, though with different out- 
comes at the nonlinear stage. The instability of the 
p-front produces smooth humps with a relatively large 
scale, about (0.1 — 0.3)mm, while the n- front appears as 
a combination of thin elongated "fingers". We focus on 
the p-front dynamics, since p-doping dominates for this 
choice of OSG sweeping more than 75% of the active ma- 
terial [m. 

A smooth shape of the p-doping front allows reducing 
the full model to a single nonlinear equation 
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where the Darrieus-Landau operator J and the Hilbert 
operator H imply products by and \k\/ik in Fourier 
space, respectively, and Xp cx Lp is the cut-off wavelength 
related to the small front thickness (detailed derivation 
will be presented elsewhere). Equation (j3j) has been 
solved numerically with initial conditions extracted from 
the experimental data at t = 10 s. The numerical model- 
ing reproduces well the characteristic shape of the p-front 
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FIG. 3. Mean front position for p- and n- fronts versus time 
obtained experimentally; predicted theoretically for the pla- 
nar fronts; and calculated numerically for the wrinkled fronts 
(see the legend). The error bars of the experimental data 
indicate difference in the positions of the fastest and slowest 
parts of the corrugated front brush. 



at later instants, e.g. at t = 70 s, see Fig. 2(d). Figs. 
2(e) and (f) show the computed electric field ahead of the 
corrugated p-type front. The observed increase in the 
electric field at the humps of the front demonstrates the 
instability mechanism, as explained earlier. The analyt- 
ical, experimental, and numerical results for the doping 
front positions are presented in Fig. 3. The analytical 
result (solid lines) shows two planar fronts accelerating 
towards each other. The acceleration takes place because 
a constant potential difference is effectively applied over 
the continuously decreasing distance between the dop- 
ing fronts [24, 25], and because the front velocities are 
controlled by the potential gradient, Eq. (1). The ana- 
lytical result relies on our experimental data [19], with 
the average hole concentration = 8.6 • lO^^m"^. 

The experiments reveal a noticeably faster propagation 
of the non-planar doping fronts compared to the analyt- 
ical result for planar fronts. In Fig. 3, open diamonds 
show the mean positions of the experimental non-planar 
doping fronts, where the "uncertainty" bars indicate the 
difference between the fastest and the slowest parts of 
the front brush. The difference can be attributed to the 
developing instability and the associated wrinkled front, 
which is not accounted for in the analytical description 
of the planar fronts. The analytical prediction for posi- 
tions of the planar fronts correlates well with the slowest 
points of the experimentally observed fronts, while the 
leading points in experiments move approximately 1.4 
times faster. In contrast, the numerical modeling does ac- 
count for the instability and shows good agreement with 
the experiments. 

The present results therefore verify that front insta- 
bility increases the doping rate. However, the instability 




FIG. 4. Experimental photographs showing the time evolu- 
tion of the p-type doping front in an LEG device with the 
initial corrugations of the anode. The corresponding time in- 
stances are indicated on the photographs. The insert shows 
characteristic shape of a strongly curved hump obtained in 
the theory. 



theory also indicates that the time and space available for 
the perturbation growth in the employed LEG devices are 
not sufficient to achieve a strongly curved front. Specifi- 
cally, the theory suggests that the instability may double 
the velocity of a curved front in comparison to a planar 
one, which is not observed in the experimental results of 
Fig. 3 discussed above. Such strong increase of the front 
velocity happens due to strongly curved humps, which 
may develop at the front (see the insert of Fig. 4). There- 
fore, to achieve a strongly curved front and a faster turn- 
on, we have triggered the instability growth by modifying 
the initial conditions so that the slowest initial stage of 
hump growth from natural "white- noise" perturbations is 
eliminated. Figure 4 presents results from an experiment 
with a corrugated pattern introduced as the surface of 
the anode (see the dotted line); the corrugation size was 
comparable to the humps observed in Fig. 2. The corru- 
gations, though negligible in comparison with the total 
(1 mm) width of the active material, were anticipated to 
give rise to the desired initial perturbations of the p-type 
front. Confirming our anticipation, we observed strongly 
curved and rapidly growing humps at the p-type front, 
see Fig. 4. The position of the humps in this experi- 
ment is also included in Fig. 3, which shows that, in 
devices with a corrugated electrode surface, the doping 
rate may increase significantly; by a factor of two in the 
present case. We have also performed experiments with 
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corrugations included into both the cathode and anode, 
and, as expected, we observed intensified front dynamics 
and much faster turn on. By introducing corrugations 
on the cathode we also find that the initial delay in the 
development of the n-type front was reduced. 

To summarize, we have demonstrated theoretically and 
experimentally that electrochemical doping fronts are un- 
stable. The obtained instability distorts the fronts and 
increases the doping rate considerably. This understand- 
ing of new fundamental properties of organic semicon- 
ductors was utilized in the design of light-emitting elec- 
trochemical cells with distinctly improved turn-on times. 
We further propose that similar design principles can be 
employed in other kinetically limited electrochemical de- 
vices, such as actuators [26[ and transistors 
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